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We report on an observation of propagating compression waves in a quasi-two-dimensional mono- 
layer of apolar granular rods fluidized by an upflow of air. The collective wave speed is an order 
of magnitude faster than the speed of the particles. This gives rise to anomalously large number 
fluctuations AAT ]s[^-^2±q.qa ^ ^^^-^^^ 

are greater than ordinary number fluctuations of N^^'^ . We 
characterize the waves by calculating the spatiotemporal power spectrum of the density. The po- 
sition of observed peaks, as a function of frequency uj and wavevector /c, yields a linear dispersion 
relationship in the long-time, long- wavelength limit and a wavespeed c = uj/k. Repeating this anal- 
ysis for systems at different densities and air speeds, we observe a linear increase in the wavespeed 
with increasing packing fraction with no dependence on the airflow. Although air-fluidized rods self- 
propel individually or in dilute collections, the parallel and perpendicular root-mean-square speeds 
of the rods indicate that they no longer self-propel when propagating waves are present. Based on 
this mutual exclusivity, we map out the phase behavior for the existence of waves vs self-propulsion 
as a function of density and fluidizing airflow. 

PACS numbers: 45.70.-n, 47.55.Lm; 87.18.Gh, 47.54.r, 05.65.b 



I. INTRODUCTION 

Macroscopic granular systems have been studied ex- 
tensively due to their ubiquity in industrial applications 
as well as the ease with which they may be visualized [l]- 
[3]. Furthermore, air-fluidized grains have demonstrated 
a remarkable ability to serve as both a system analogous 
to Brownian particles ^ as well as a proxy for colloidal 
and molecular glasses [51 |6]. Idealized granular media 
consist of spherical particles whereas real systems will 
have both polydispersity and anisotropic particle shape. 
The role of particle shape is now a topic of increasing at- 
tention. Recent experiments on vibrated rods [7HTT] and 
rods in a hopper [12 have demonstrated the great vari- 
ety of phenomena - swirling, vortices, pattern formation 
- that arise due to particle anisotropy alone. 

Individual air- or vibro-fluidized granular rods have 
also been observed to self-propel, in the sense of mov- 
ing preferentially along their long axis p!3HT5] . The idea 
of active self-propulsion has garnered great theoretical 
interest [161 III] due to the potential of unifying a broad 
spectrum of systems over a huge range of length scales: 
from biological systems like bacteria [18H21], fish \22\ . 
locusts [23] , birds [24l [25] , to physical systems such as 
agitated granular materials [Tj |TT] and human or ani- 
mal traffic [26l[27]. In this context, self-propelling at a 
microscopic scale gives rise to dramatic collective behav- 
ior at the macroscopic scale including the emergence of 
a dynamic broken-symmetry state in which all the par- 
ticles move in the same spontaneously-chosen direction 
as a coherent flock [28-30 . These systems are also pre- 
dicted to exhibit collective motion, propagating waves, 
and anomalously large number fluctuations [3TH35] . De- 
spite the wealth of theory and simulation, very few exper- 
iments have been conducted due to both the difficulty of 
tracking and analyzing biological systems as well as find- 
ing purely physical systems in which particles self-propel. 



Of note, giant number fluctuations were observed in a 
vertically-vibrated monolayer of granular rods possess- 
ing nematic order [T3|. Additionally, collective behavior 
in the form of whorls, jets, and vortices - in which the 
collective motion was an order of magnitude faster than 
individual particle motion - has been observed in dense 
populations of bacteria [19 . 

For this paper, we extend our earlier experiments on 
individual and dilute systems of air-fluidized rods [15] in 
order to search for giant number fluctuations and col- 
lective flock-like motion. In spite of considerable effort, 
we never found conditions under which flocking occurred. 
We do observe giant number fluctuations, however, and 
we additionally find that they appear to be coupled to 
the unexpected existence of density waves that propa- 
gate ballistically across the system. To begin, in Sec- 
tion [nj we introduce the apparatus and system under 
study. The paper will first examine a single system at 
fixed density of 64% packing fraction. We describe our 
observation of propagating compression waves in Sec- 



tion III The wavespeed is an order of magnitude larger 
than the root-mean-square speed of the particles; this 
separation of timescales results in giant fluctuations in 
the local number density, discussed in Section [TV| Lastly, 
in Section [V| we quantify the waves by calculating the 
spatiotemporal power spectrum of the density and ex- 
tracting a dispersion relationship for the waves. We ex- 
tend the experiments by examining how the wavespeed 
varies with density and airflow and mapping out phase 
behavior. 



II. EXPERIMENTAL DETAILS 

We study a monolayer of right cylindrical acetal dowel 
rods - length 1.27 cm, diameter 0.24 cm, aspect ratio 5, 
and mass 0.076 g - fluidized by an upflow of air. Except 
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FIG. 1: (Color online) A propagating wave in a monolayer of gas-fluidized rods. The system has diameter 28.6 cm. The rods 
are 1.27 cm long, have an aspect ratio of five, occupy 64% projected-area packing fraction, and are fluidized by an upflow of 
air at 220 cm/s. The time between images, moving left to right, is 0.15 s. The red arc is moved at constant velocity 28 cm/s 
as a guide to the eye. 



when otherwise noted, we analyze a system of 1353 rods 
occupying packing fraction 64%, fluidized by a spatially 
and temporally uniform upflow of air at speed 220 cm/s 
as measured by a hot-wire anemometer. The fluidizing 
air speed is low enough that rods do not overlap out-of- 
plane and the system remains a monolayer. 

The apparatus and fluidization method is identical to 
that of Ref. [15]. The apparatus is a rectangular wind- 
box, 0.5 X 0.5 X 1.2 m^, positioned upright. A circular 
brass testing sieve with mesh size 150 /im and diame- 
ter 30.5 cm rests horizontally on top. In order to pre- 
vent particles from becoming trapped in a small groove 
around the edge of the sieve, we place a 1 cm norprene 
tube around its inside edge. A blower attached to the 
windbox base provides vertical airflow perpendicular to 
the sieve. Raw video data of the fluidized particles is cap- 
tured for 10 minutes at 120 frames per second by a digital 
camera mounted above the apparatus. Post-processing of 
the video data is accomplished in Lab VIEW, using the 
same tracking programs used in Ref. [15j when appropri- 
ate. For the analysis of the waves in Section [V| we treat 
the video as a binary density map t) rather than 

individually track particles. 

In order to emphasize the role that shape plays here, 
we analyze a companion system of 400 bidisperse plastic 
spheres - diameters 0.64 cm and 0.95 cm - chosen to 
have the same packing fraction 64% as the rod system. 
The fluidization speed, 260 cm/s, for the sphere system 
is chosen so that the the spheres and rods have roughly 
the same rms speed. This is below the terminal falling 
speed, so the spheres roll without slipping or levitating. 

Prior to focussing on the 1.27 cm rods with aspect 
ratio 5, we also explored a number of other granular sys- 
tems. This includes apolar rods with aspect ratios of 4, 
9, 14, and 30; polar rods having a distinct head and tail 
with aspect ratio 2.5 and 14; a range of airspeeds from 
10 to 500 cm/s; and densities from 10% to 85% area 
packing fraction. We also examined such systems in a 
circular race track with width that varied from approxi- 
mately three to ten particle widths, by placing a cylinder 



concentrically onto the sieve. In no case did we observe 
flocking. 



III. PROPAGATING WAVES 

We observe a propagating compression wave instability 
for dense collections of fluidized rods. This phenomenon 
is unique to rods; we do not observe compression waves 
for spheres. A time series of images, each separated by 
0.09 s, shown in Figure [l] depicts one such propagating 
compression wave for a system at 67% packing fraction. 
The red arc moves at constant velocity 28 cm/s with the 
wavefront and serves as a guide to the eye. When parti- 
cles are compressed together, light is specularly reflected 
from the interface and the wavefront appears brighter in 
the images and the less dense rarefaction zones darker. 
An example video clip is available on-line [36^ . 

As can be seen in videos, e.g. [36], the waves travel 
much faster than the particles themselves. The waves 
propagate through the medium regardless of the local 
ordering with some preference for moving along the local 
splay direction. To get an estimate for the wavespeed 
without any detailed processing, we look at spacetime 
plots of a single line of pixels extracted from the video 
data, as shown in Fig. [2] The waves are visible on 
the spacetime plot as diagonal lines. From the slopes 
of these lines, we obtain a range of wavespeeds from 
10 to 30 cm/s, in agreement with the moving arc in 
Fig. [l] From position- versus-time data for the individual 
rods, we calculate the root-mean-square particle speed 
of ~1 cm/s. Thus, the waves move an order of magni- 
tude faster than the particles themselves and an order 
of magnitude slower than the fluidizing airflow. This is 
similar to what has been observed for collective motion 
of bacteria in which the collective behavior is an order of 
magnitude faster than the individual particle speeds [19] . 
However, both a hydrodynamic theory of particles ad- 
sorbed on a substrate [Sj as well as a microscopic the- 
ory of self-propelling rods [35] predict that propagating 



FIG. 2: Spacetime plot of a line of pixels taken from a video 
of air-fluidized rods as specified in Section [ill 



waves in those systems have a speed c equal to Vrms- 

Predictions of waves for self-propelling rods typically 
require that long-range order or a dynamic broken sym- 
metry state be present in the system. By observation, 
neither appears to develop for our systems. The emer- 
gence is suppressed by the propagation of the compres- 
sion waves which disrupt the local ordering and any co- 
herent motion developing between neighboring particles. 
However, although there is no long-range order, localized 
domains of nematic ordering are observed and the parti- 
cles in these regions tend to move in the same direction 
in a bulk sense - we observe localized correlated motion. 
To quantify the extent of these local domains, we calcu- 
late the ensemble- and time-averaged orientational order 
parameter (|e*^^|) for n = 2, 3, 4. For a fixed packing 
fraction and airflow, we vary the neighborhood size over 
the which the ensemble average of orientations is taken. 
The neighborhood is a square subregion with side lengths 
/ X The result for a system at 64% packing fraction is 
shown in Fig. [3] where the subregion side length has been 
scaled by the rod length. If a system possesses an n- 
fold symmetry within the subregion, the quantity will be 
equal to 1. The results show that there is substantial 
nematic (n=2) local ordering with values of the order 
parameter ranging from 0.6 to 0.7 in a neighborhood of 
two particle lengths, indicating the tendency for neigh- 
boring rods to align along their long axes; the values for 
3- and 4-fold symmetry are also appreciable. The quan- 
tity decays to approximately 0.1 when we average over 
the entire system. A fit to an exponential yields a 1/e 
correlation length of approximately 8 particle lengths for 
nematic order - nearly half the system size - indicating 
medium-range order, intermediate between particle size 
and system size. 

We have no firm explanation for the existence of prop- 
agating waves that are an order of magnitude faster than 
the particle speed. One possibility is based on a Peclet 
number [19], defined as the Peclet number is the ra- 
tio of advection to diffusion UL/D - [/, the particle 
speed; L, the particle length; and the diffusion con- 
stant found from the linear regime of the particle mean 



FIG. 3: (Color online) Time- and ensemble- averaged orien- 
tational order parameter (|e*^^|) for n= 2, 3, 4 for a system 
of air-fluidized rods at 64% packing fraction. The ensemble 
average is taken over all particle orientations within a square 
subregion with side lengths / x /. The order parameter is 
plotted against this side length scaled by the rod length. 



square displacements obtained by individually tracking 
the rods [15 . For this argument, as an individual rod 
moves preferentially parallel to its long axis, the sur- 
rounding granular fluid becomes entrained and gives rise 
to advection flows. The entrained fluid causes neighbor- 
ing rods to prefer to align in the direction of the advective 
flow and travel in the same direction. This is one way to 
understand the localized correlated motion observed in 
our system. When these localized domains collide with 
one another, the particles are constrained from moving 
but the entrained fluid flows past the particles, pushing 
the neighboring domain and resulting in a propagating 
wave. For our system, the Peclet number UL/D - using 
U ^ 1 cm/s and D = 0.128 cm^/s - is approximately 10, 
suggesting that advection could play a role. A second 
possible explanation would entail an inertial argument 
in which the rods attempt to self-propel along their long 
axis but are constrained by the reduced free area avail- 
able to them. Thus, during a collision with another rod, 
there is an additional transfer of momentum that pushes 
the second rod into a third rod, giving rise to a cascade 
of collisions that result in a propagating wave. 



IV. NUMBER FLUCTUATIONS 

One interesting consequence of the propagating waves 
moving much faster than the individual rods is that com- 
pressed and dilute regions coexist for periods of time long 
compared to the time it takes the waves to propagate. 
Such persistent fluctuations in the density suggest the 
presence of anomalously large fluctuations in the local 
number density for these systems, which we now calcu- 
late. For a thermal system, the quantity AN/ ^J (N), 
where {N) is the mean number of particles in a subregion 
and AN is the standard deviation, should be a constant. 
Thus, thermal behavior will be characterized by a hori- 
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FIG. 4: (Color online) Magnitude of number fluctuations nor- 
malized by v^(AO versus average number of particles in a sub- 
region of the system. The circles (green) are for a bidisperse 
collection of spheres. The squares (blue) are for a collection 
of rods. Both systems occupy 64% areal packing fraction and 
are fluidized at 220 cm/s. The vertical dashed lines corre- 
spond to the total number of particles in each system: 1353 
rods and 400 spheres. The dashed line is a power-law fit, 
which corresponds to AN - (iV)0-72±o.o4^ 



zontal line on a plot of AN/ ^J {N) versus (A^), whereas 
so-called giant number fluctuations will be characterized 
by a non-zero slope and a magnitude larger than one. 

To quantify number fluctuations, we take a 10-minute 
video of the system at a given packing fraction and air- 
flow. We select a square subregion of interest with side 
length / and count the number of particles within that 
region for each frame. From this time series, we calculate 
(A^) and AA^. We then repeat this procedure for system- 
atically larger regions of interest, ranging from a single 
pixel up to roughly half of the system size. 

Results for a system at 64% packing fraction, fluidized 
at 220 cm/s, are shown in Figure Ul plotting AN/ (N) 
versus (N) for each subregion. The number fluctua- 
tions for bidisperse spheres show ordinary behavior, with 
AN/ y^J^ ^ 1.5 being constant over the entire range 
of subregion sizes. Rods, however, show number fluctua- 
tions with both an exponent and magnitude larger than 
thermal expectations, indicating giant fluctuations in lo- 
cal number density. The large (N) fall-off is a finite size 
effect; the vertical dashed lines in Fig. |4] indicate the to- 
tal number of particles in each system. At its maximum 
prior to this fall-off, the fluctuations in local number den- 
sity of the rods are larger than those for spheres by more 
than a factor of 3. By fitting a power law to the rod 
data, we find that number fluctuations for the rods scale 
as AN - (7V)0-72±o.04^ 

Although the maximum in the number fluctuations is 
due to finite size effects, it serves as a useful benchmark 
to quantify the spatial extent of local number fluctua- 
tions. We accomplish this by converting the value of AA^ 
into an effective range of packing fractions. The value 
of A^ + AA^ at the maximum corresponds to a range of 
packing fractions (j) ~ 64% ± 20%. Thus, as the wave 
propagates through the system, particles are compressed 
by the front up to 84% while the rarefaction zones are 



diluted to approximately 44%. 

Figure[4]is strikingly similar to the number fluctuations 
obtained for rods fluidized by vertical vibrations [13 . In 
that experiment though, the large voids responsible for 
giant number fluctuations persisted for very long times, 
indicated by a logarithmic decay in the local density au- 
tocorrelation function. In our system, the density au- 
tocorrelation function decays to oscillations about zero 
within a few seconds. These two systems exhibit qual- 
itatively similar behavior - giant number fluctuations - 
caused by different mechanisms. The common factor be- 
tween the two systems is that the rods self-propel at low 
densities. 



V. WAVES AND DENSITY 

Propagating waves have been predicted for collections 
of self-propelling particles interacting via volume exclu- 
sion with no emergent broken-symmetry state [35 , as 
well as more generally for particles adsorbed on a surface 
in a fluid background i37|. In order to facilitate a the- 
oretical description of our observations and understand 
the relative roles of particle shape and self-propelling, 
we must know how wavespeed and particle speed change 
with packing fraction and airflow. We do this by cal- 



culating a dynamic structure factor in subsection |V A 



We repeat the analysis for systems at different packing 



fractions and airflows in subsection |VB[ Lastly, in sub- 
section V C[ we map out the phase behavior of the rods. 



A. Density Power Spectrum 

In this subsection, we obtain more detailed quanti- 
tative information about the waves by calculating the 
spatiotemporal power spectrum C{uj^\k\) of the density 
p{x^y^t) - where uj is frequency and \k\ is wave-vector. 
Rather than compute the actual density, we simply use 
the 0-255 grayscale pixel values directly from the video 
data. As seen in Fig. [l] this may be taken as a proxy 
for the local density. Using LabVIEW's Vision pack- 
age, we first obtain the spatial Fourier transform of each 
frame in the video. Because there is no long-range order, 
the spatial Fourier transform of the rods data shows two 
rings at 4.96 cm~^ and 26.2 cm~^, corresponding 

to the long and short dimensions of the rod. The annuli 
are isotropic with respect to the polar angle indicating 
that there is no long-range order in the system. We next 
extract the time-traces for all pixel values within an an- 
nulus at fixed \k\. We then calculate the temporal power 
spectrum for each pixel time-trace at fixed \k\ and av- 
erage over all pixels to obtain the spatiotemporal power 
spectrum C{uj^\k\). The magnitude of the power spec- 
trum as both u and \k\ is simply the sum of the 
grayscale values of each frame in the video averaged over 
time. We normalize the power spectrum magnitude by 
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FIG. 5: (Color online) Spatiotemporal power spectrum of the 
density C{uj,k) for a) a collection of rods and b) a collection 
of bidisperse spheres, both occupying 64% packing fraction. 
Each curve corresponds to a different wave-vector value as 
shown by the color legend in a). The DC-component for each 
curve is plotted to the left of the axis break. 



this quantity so that the rods and spheres data can be 
more directly compared. 

Slices of C{uj^\k\) for fixed \k\ are plotted versus uj 
in Fig. (sja) for rods and jsjb) for spheres, both at 
64% projected area fraction. The range of wave- vectors 
shown correspond to one-fourth the system size to ap- 
proximately one particle length. The salient feature of 
C{uj,\k\) in Fig. (sja) is a single peak at a given in- 
dicating a traveling excitation. For larger wave-vectors 
outside of the hydrodynamic limit, the peak is no longer 
observed. In contrast, the same analysis for spheres in 
Fig. [sjb) shows no such peaks. 

By plotting the location of peak frequency versus 
we can construct the dispersion relationship for the 
waves, as shown by the symbols in Fig. [6j On this 
same figure, we also superpose a colorized contour plot of 
C{uj^\k\). In the limit of small wavevector and frequency, 
the dispersion relation is linear and we may therefore ex- 
tract a wave speed c = uj/\k\ from the slope. For this par- 
ticular example, c = 20.8 cm/s. This is consistent with 
the video data which show the waves propagating across 
the system in 1-2 seconds as well as the value obtained 
from spacetime plots. Again, we note that the wave speed 
is an order of magnitude larger than the rms speed of 
the particles, Vrms ^ 1 cm/s, and an order of magnitude 
smaller than the fluidizing airflow ~ 200 cm/s. 




FIG. 6: (Color online) Contour plot of the spatiotemporal 
power spectrum of the density for rods occupying 64% pack- 
ing fraction and fluidized at 220 cm/s. The power spectrum 
amplitude is given in the color legend. The ocations of the 
peaks (circles) form the dispersion relationship for the propa- 
gating waves. The slope of the linear region is the wave speed 
c 20.8 cm/s. 



B. Wavespeed 

Now we may measure the wavespeed as a function of 
both packing fraction and airspeed, using the spatiotem- 
poral power spectra as illustrated in the examples of 
Figs. [5][6j For comparison, we also use particle tracking 
to measure the rms speeds parallel and perpendicular to 
the long axis of the rods. The results are displayed in 
the four plots of Fig. |8] At fixed air speed, the wave 
speed increases slightly as we increase the packing frac- 
tion of the system as seen in Fig. [7|^a). The increase ap- 
pears to be roughly linear, although the dynamic range 
of our data is limited. This linear trend is seen for sev- 
eral other fixed airflows that we analyzed. By contrast, 
if we fix the packing fraction and increase the fluidizing 
airflow, the wavespeed shows no dependence on airflow, 
Fig.[7|;b). From Fig. [7|;c), we see that the particle speeds 
increase slightly with increasing density; from (d), we see 
a similar increase in particle speed with increasing airflow 
except at the highest airflows. For high airflows, out-of- 
plane motion dominates and in-plane motion slows down. 
For most densities, the parallel self-propelling speed is 
slightly larger than the transverse speed, although the ef- 
fect is smaller than that observed for dilute systems [15] . 
Although the rods are self-propelling, we see that the 
side-to-side motion induced by the propagating waves 
causes v:^^^ to be nearly equal to Vrms- 

For all densities and airspeeds analyzed, the wavespeed 
remains an order of magnitude larger than the parti- 
cle speed and the components of the particle speed re- 
main nearly equal. This indicates that when propagat- 



6 




0.55 0.6 0.65 200 
Packing Fraction 



250 300 
Airflow [cm/s] 



FIG. 7: (Color online) Wavespeed as a function of a) packing 
fraction and b) airflow and root-mean-square parallel (open 
symbols) and perpendicular (closed symbols) particle speeds 
as a function of c) packing fraction and d) airflow. 
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FIG. 8: (Color online) Wavespeed (blue circles), c — uj/k^ and 
root-mean-square parallel (red open diamonds) and perpen- 
dicular (black closed diamonds) particle speeds as a function 
of packing fraction for rods fluidized at 280 cm/s. The shaded 
green region is where out-of-plane motion becomes substan- 
tial. 



ing waves are present, the particles no longer self-propel. 
This is demonstrated more powerfully in Fig. |8) where 
we extend the particle rms speed measurements to lower 
packing fractions. For this particular density and airflow, 
waves are no longer observed below 50% packing frac- 
tion. And below this packing fraction, we see that the 
parallel and perpendicular particle speeds become appre- 
ciably different. Thus, we conclude that self-propulsion 
and waves are exclusive of one another. 



FIG. 9: (Color online) Behavior of fluidized rods as a function 
of fluidizing airflow and packing fraction. The diagram in a) 
is for a large system with diameter 29 cm; b) is for a smaller 
system with diameter 14 cm. The light blue cross in a) are 
the conditions for the analysis of Sections [V|\. The orange 
cross in b) represents the conditions for the experiment of 
Ref. Jj]. The vertical dashed line in both plots corresponds 
to the density above which the rods are uniformly distributed 
about the system. 



C. Phase Behavior 

Since we observed that self-propelling and waves are 
mutually exclusive, we want to determine the conditions 
for which a system of rods exhibits each phenomenon. 
Here, we qualitatively map out the phase behavior of flu- 
idized rods as a function of packing fraction and fluidizing 
airflow. Fig. |9] We explored a range of densities from a 
single particle up to 75% packing fraction and fluidiza- 
tion speeds from 150 cm/s to 500 cm/s. We first place the 
desired number of rods into the system with no air flow. 
We then turn the airspeed to 150 cm/s. We increase the 
airspeed slowly by increments of roughly 10 cm/s. After 
each increase in the airflow, we wait 1 minute to ensure 
the system has reached a steady state and then character- 
ize the behavior observed before increasing airflow once 
more. There are strong particle- wall interactions that 
cause the rods to cluster in the center of the system at 
low densities. Thus, the exact shape of the phase bound- 
aries system-size-dependent. As such, we show a phase 
diagram for the full system size - diameter 28.6 cm - in 
Fig. |9ja) as well as the system size used in Ref. - di- 
ameter 17 cm - in Fig. [9]^b). Because packing fraction is 
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ill-defined unless the particles are uniformly distributed 
across the system, we denote the packing fraction above 
which the rods are uniformly dense as a vertical dashed 
line in both Figs.|9|a) and (b). 

For all packing fractions, there is a threshold airflow 
below which the rods do not have enough energy to move. 
This is the lower boundary of the phase diagrams shown 
in Fig. [9j Conversely, above some airflow, the rods gain 
enough energy to lose contact with the substrate and 
overlap one another out-of-plane. This marks the upper 
boundary of the phase diagrams which is determined by 
when we first observe two rods overlapping. Within these 
boundaries, we observe two distinct behaviors depending 
on both packing fraction and airflow. For dilute systems, 
the rods smoothly self-propel with very little out-of-plane 
motion. The experimental conditions for our previous 
work are shown as the solid orange cross in Fig. [oj^b) . 
As density is increased, there emerges an instability for 
compression waves to propagate. The onset of waves is 
chosen to be when ripples first appear to the eye. 

VI. CONCLUSION 

We have investigated the behavior of dense collections 
of rods fluidized by an upflow of air. For a large range 
of densities and airflows, the monolayer of rods becomes 
unstable to the propagation of compression waves. The 
waves travel at a speed an order of magnitude larger 
than the root-mean-square speed of the particles them- 
selves. As the compressed and rarefied regions relax after 



a wave's passage, the separation in timescales gives rise 
to anomalously large fluctuations in the local number 
density. Unlike giant number fluctuations expected for 
swarms and seen in vibrated-bed experiments, number 
fluctuations in our system are short-lived. Nevertheless, 
the particles in both experiments are self-propelling rods 
and these experiments suggest that self-propelling can 
be used as a useful concept to unify both physical and 
biological systems. 

To ascertain the roles that self-propulsion and shape 
play in the collective wave behavior we observe, we quan- 
tified the waves by calculating the dynamic structure 
factor for the rods. From the position of peaks in this 
function, we extracted a dispersion relationship and, in 
the hydrodynamic limit, the wavespeed. Repeating our 
analysis for systems of increasing density and constant 
airflow, we found that the wavespeed is roughly linear 
with particle density and independent of airflow. Fur- 
thermore, for packing fractions at which waves are ob- 
served, the parallel and perpendicular rms particle speeds 
are the same, indicating that waves and self-propelling 
behavior are exclusive of one another for our system. 
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